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Magnetic  Field  Interactions 
with  Material  Discontinuities 


by 

W.  Lord 


Abstract 

All  methods  of  nondestructive  testing  (NDT)  rely  for  their  operation 
on  the  interaction  of  an  energy  source  with  defects  in  the  material  under 
test.  Such  interactions  are  extremely  difficult  to  model  analytically 
because  of  the  awkward  boundaries  associated  with  realistic  defect  shapes 
and  the  inherent  nonlinear  behavior  of  many  material  properties.  The  major 
contribution  of  this  project  work  has  been  to  show  that  finite  element 
analysis  techniques  can  be  used  to  model  electromagnetic  energy/defect 
interactions  associated  with  active  leakage  field,  residual  leakage  field, 
eddy  current  and  pulsed  eddy  current  NDT  methods. 

Developed  code  has  been  used  to  study  probe  design,  to  simulate 
testing  geometries  too  difficult  to  replicate  in  a  laboratory  environment 
and  to  aid  in  the  development  of  defect  characterization  algorithms. 
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Electromagnetic  methods  of  nondestructive  testing  are  used  widely 

throughout  the  aerospace,  transportation,  ordnance,  metals  and  energy 

industries  for  detecting  defects  in  critical  metal  parts.  A  major 

stumbling  block  in  the  development  of  these  techniques,  particularly  with 

regard  to  the  study  of  automated  defect  characterization  schemes  and 

associated  imaging  systems,  has  been  the  lack  of  a  suitable  model  capable 

of  predicting  the  inherently  complex  field/defect  interactions. 

Over  the  past  decade,  with  Army  Research  Office  support,  the  NDT 

laboratory  at  Colorado  State  University  has  pioneered  in  the  use  of  finite 

element  analysis  techniques  to  solve  the  modeling  problem  for  active  and 

residual  leakage  field  as  well  as  eddy  current  NDT  methods^. 

With  the  advent  of  the  digital  computer,  increasing  emphasis  has  been 

placed  on  the  numerical  solution  of  partial  differential  equations.  In  the 

2 

late  1960's  Erdelyi  showed  how  the  finite  difference  method  of 

approximating  partial  derivatives  could  be  applied  to  the  prediction  of 

electromagnetic  fields  in  electrical  machinery.  At  about  the  same  time 
3  4 

Winslow  and  Silvester  and  Chari  demonstrated  that  finite  element 
techniques,  based  on  variational  calculus,  could  be  applied  to  the  same 
problems.  A  period  of  acrimonious  debate  followed  in  the  power  apparatus 
and  systems  literature  as  to  the  relative  merits  of  the  two  approaches^. 

Based  on  the  author's  experience  in  the  early  197U's  of  applying 
finite  element  analysis  to  the  study  o  electromagnetic  fields  in  magnetic 
structures'*^,  it  became  clear  that  such  techniques  would  be  ideally  suited 
to  the  study  of  electromagnetic  NDT  phenomena  because  of  the  ease  with 
which  the  awkward  defect  boundaries  could  be  handled. 

With  a  small  "seed-money"  grant  from  the  Colorado  Energy  Research 


Institute  in  1974,  initial  work  was  undertaken  to  develop  finite  element 
code  for  the  study  of  leakage  fields  around  defects  in  a  cylindrical  shaped 
bar  carrying  an  axial  dc  excitation  current  -  the  "active”  case.  Based  on 

Q 

promising  preliminary  results  from  the  CERI  study  ,  a  proposal  was 

submitted  to  the  Army  Research  Office  in  1975.  The  major  objectives  of  the 

ARO  study  were  to  a)  extend  the  preliminary  results  to  a  wider  variety  of 

defect  shapes,  b)  determine  the  feasibility  of  extending  the  finite 

element  code  to  the  residual  leakage  field  (magnetic  particle)  NDT 

situation  and  c)  explore  the  possibility  of  applying  finite  element 

analysis  techniques  to  a  study  of  eddy  current  NDT  phenosiena .  Following 

9-14 

successful  completion  of  these  objectives  ,  further  work  was  initiated 
via  an  ARO  follow-on  proposal  aimed  at  extending  the  studies  to: 

a)  The  examination  of  an  unusual  leakage  field  reversal  phenomena^^, 

b)  three  dimensional  magnetostatic  finite  element  code^^, 

17  18 

c)  pulsed  eddy  current  phenomena  '  ,  and 

d)  deconvolution  effects  ^°X9  element  measurements  of 

magnetostatic  leakage  fields^  . 

In  short,  Army  Research  Office  support  of  studies  in  the  Applied 

Magnetics  and  NDT  Laboratory  at  Colorado  State  University  relating  to 

electromagnetic  field/defect  interactions  have  clearly  shown  the 

feasibility  and  usefulness  of  applying  finite  element  analysis  techniques 
20  21 

as  a  modeling  tool  *  .  In  addition,  funding  from  the  Electric  Power 

Research  Institute  has  demonstrated  the  applicability  of  such  a  modeling 

22-30 

tool  to  practical  NDT  problems  of  interest  to  the  nuclear  industry 

All  this  work  has  clearly  shown  the  close  interrelationship  between 

models  of  active,  residual  and  eddy  current  NDT  phenomena,  thus  paving  the 

way  for  a  unified  approach  to  the  problem  of  modeling  electromagnetic 

21 

field/defect  interactions  . 


In  carrying  out  the  modeling  studies  it  has  also  become  apparent  that 

the  finite  element  code  has  a  number  of  potential  uses: 

a)  An  an  "experimental  model"  for  simulation  of  electromagnetic  NDT 
situations  too  difficult  or  expensive  to  replicate  in  a  laboratory 
environment . 

b)  As  an  aid  to  the  physical  understanding  of  the  interactions 
between  electromagnetic  fields  and  defects  in  the  part  under  test. 

c)  As  a  design  tool  fox  the  study  of  alternative  probe  geometries. 

d)  As  a  training  mechanism  for  the  development  of  automated  defect 
characterization  schemes. 

e)  As  a  technique  for  estimating  the  bulk  properties  (conductivity 
and  permeability)  of  materials  under  test.  The  procedure  consists 
of  iterating  the  finite  element  code  with  different  property 
values  until  the  code  predictions  agree  with  experimental 
observations . 

2.  SiiwiiiaT-Y  of  Results 

Finite  element  code  has  been  developed  for  the  following 

electromagnetic  NDT  testing  situations: 

a)  Active  leakage  fields: 

1.  2-dimen8ional . 

2.  Axisymmetric . 

3 .  3-dimens iona 1 . 

b)  Residual  leakage  fields: 

1.  2-dimensional. 

c)  Single  frequency  eddy  current: 

1 .  2-dimensional . 

2.  Axisymmetric. 

3.  3-dimensional. 

d)  Pulsed  eddy  current: 

1.  Axisymmetric. 

In  addition,  experimental  work  has  been  carried  out  using  the  test  rig 
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shown  schematically  in  Figure  1.  This  has  allowed  confirmation  of  the 

developed  finite  element  code  and  led  to  some  novel  studies  of  leakage 

15  19 

field  reversal  phenomena  ,  Hall-element  deconvolution  and  the  imaging  of 

31 

magnetostatic  leakage  fields  and  eddy  current  signals  (see  Figs.  2  and 

3).  Initial  work  has  also  begun  on  the  application  of  finite  element  code 

32 

to  the  simulation  of  ultrasonic  NDT  phenomena 
3 .  Publications 


From  April  1980  to  September  1984  the  following  journal  publications 

have  been  written  with  support  from  ARO: 

A  survey  of  electromagnetic  methods  of  nondestructive  testing.  Chapter 
3  in  the  text  Mechanics  of  Nondestructive  Testing .  edited  by  W. 
W.Stinchcomb,  Plenum  Press,  1980,  pp.  77-100. 

Finite  element  analysis  of  eddy  current  phenomena.  Materials 
Evaluation.  Vol.  38,  No.  10,  Oct.  1980,  pp.  39-43. 

Development  of  theoretical  models  for  NDT  eddy  current  phenomena,  in 
Eddy  Current  Characterization  of  Materials  and  Structures.  ASTM  STP 
722,  G.  B.  Bimbaum  and  G.  Free,  Eds.,  ASTM,  1981,  pp.  5-21. 

Developments  in  the  finite  element  modeling  of  eddy  current  NDT 
Phenomena  at  CSU,  ibid,  pp.  357-361. 

Numerical  modeling  of  electromagnetic  NDT  phenomena,  in  New  Procedures 
in  Nondestructive  Testing.  P.  Holler,  Editor,  Springer-Verlag ,  1983, 
pp.  461-470. 

3-D  finite  element  prediction  of  magnetostatic  leakage  fields,  IEEE 
Transactions  on  Magnetics.  Vol.  MAG-19,  No.  5,  September  1983,  pp. 
2260-2265. 

Eddy  current  probe  design  using  finite  element  analysis.  Materials 
Eva lust ion.  Vol.  41,  No.  12,  November  1983,  pp.  1389-1394  (with  N. 
Ida  and  R.  Palanisamy) .  ASNT  1984  Achievement  Award. 

Applications  of  numerical  field  modeling  to  electromagnetic  methods  of 
nondestructive  testing,  IEEE  Transactions  on  Magnetics.  Vol.  MAG-19, 
No.  6,  November  1983,  pp.  2437-2442. 

Solution  of  linear  equations  for  small  computer  systems.  International 
Journal  for  Numerical  Methods  in  Engineering.  Vol.  20,  No.  4,  April 
1984,  pp.  625-641  (with  N.  Ida). 

Superposition  of  eddy  current  probe  signals.  Materials  Evaluation, 
Vol.  42,  No.  7,  June  1984,  pp.  930-933  (with  D.Home  and  S.  Udpa). 


Finite  element  modeling  ot  pulsed  eddy  current  phenomena,  in  Review  of 
Progress  in  Quantitative  Nondestructive  Evaluation.  D,  0.  Thompson  and 
D.  E.  Chimenti,  Eds.,  Plenum,  1984,  pp.  561-568  (with  B.  Allen). 

Deconvolution  of  defect  leakage  field  profiles  obtained  by  using  Hall 
element  probes,  ibid,  pp.  855-862  (with  L.  Srinivasan) . 

In  addition,  the  following  conference  presentations  have  been  made: 

Modeling  residual  and  active  leakage  fields,  an  invited  paper 
presented  at  the  ASMT  Fall  Conference.  Houston,  Oct.,  1980. 

Recent  developments  in  electromagnetic  NDT  methods,  an  invited  paper 
presented  at  the  Symposium  on  Novel  NDE  Methods  for  Materials.  AIME 
Annual  Meeting,  Dallas,  Texas,  February  1982. 

NDE  Education/Training  for  Engineers,  ASME-PVP  Annual  Meeting. 
Orlando,  June  1982. 

Numerical  modeling  of  electroisagnetic  NDT  phenomena,  an  invited  paper 
presented  at  the  German-Aaerican  Workshop  on  NDT  Research  and 
Deve lopment .  Saarbrucken,  V.  Germany,  August  1982. 

Eddy  current  techniques  and  their  potential  for  determining  material 
properties,  an  invited  paper  presented  at  the  Symposium  on 
Nondestructive  Methods  for  Material  Property  Determination.  Hershey, 
PA,  April  1983. 

Applications  of  numerical  field  modeling  to  electromagnetic  methods  of 
NDT,  an  invited  paper  presented  at  COMPUMAG .  Genoa,  Italy,  May  1983. 

Residual  and  active  magnetic  leakage  field  modeling,  presented  at 
QualTest  II.  Dallas,  October  1983. 

Numerical  modeling  of  eddy  current  NDE,  an  invited  paper  presented  at 
the  29th  Annual  Conference  on  Magnetism  and  Magnetic  Materials. 
Pittsburgh,  November  1983. 

The  case  for  numerically  modeling  electromagnetic  NDT  phenomena, 
presented  at  the  Seminar  on  Nondestructive  Inspection  of  Ferromagnetic 
Materials .  Dresser  Industries,  Houston,  March,  1984. 

A  new  technique  for  modeling  hysteresis  loop  phenomena,  presented  at 
an  Eddy  Current  Seminar,  Rutherford  Appelton  Laboratory,  Abingdon, 
Oxford,  U.K.,  April  1984. 

NDT  Education  at  Colorado  State  University,  The  case  for  numerical 
modeling  of  electromagnetic  NDT  phenomena,  both  papers  presented  at 
the  ASNT  Spring  Conference.  Denver,  May  1984. 


Nondestructive  evaluation  in  electrical  engineering,  presented  at  the 


Diffusion,  waves,  phase  and  eddy  current  imaging,  presented  at  ttie 
Review  of  Progress  in  Quantitative  NDE.  La  Jolla,  CA,  July  1984  (with 
L.  Udpa) . 

A  finite  element  formulation  for  ultrasonic  NDT  modeling,  presented  at 
the  Review  of  Progress  in  Quantitative  NDE.  La  Jolla,  CA,  July  1984 
(with  R.  Ludwig) . 

4.  Participating  Scientific  Persons^ 


Several  graduate  students  have  been  associated  with  the  project  during 
the  1980-84  time  frame  including: 


R.  Palanisamy 


S.  Udpa  (Satish) 


N.  Ida 


-  Ph.D.,  1980.  Currently  a  senior  NDT  research 
engineer  with  the  Timken  Co.,  Canton,  Ohio. 

-  Ph.D.,  1983.  Currently  an  Assistant  Professor 

of  Electrical  Engineering,  Colorado  State  University, 

-  Ph.D.,  1983.  Currently  an  Assistant  Professor 
of  Electrical  Engineering,  University  of  Akron. 


L.  Udpa  (Srinivaaan)  -  Ph.D.  student  at  Colorado  State  University. 


R.  Ludwig 


B.  L.  Allen 


S.  Heath 


C.  Wang 


-  Ph.D.  student  at  Colorado  State  University. 

-  M.S.,  1983.  Research  engineer  with  Rockwell 
International . 

-  M.S.,  1983,  Engineer  with  Hewlett  Packard. 

-  Currently  completing  his  Ph.D.  degree  at  North 
Carolina . 


V.  Lord  has  been  principal  investigator  throughout. 

5.  Future  Work 

This  project  has  clearly  shown  the  feasibility  of  applying  finite 
element  analysis  to  the  modeling  of  all  electromagnetic  NDT  phenomena. 
Additional  funding  has  been  requested  from  ARO  to  extend  this  work  to  the 
application  of  imaging  for  the  display  of  results  from  electromagnetic  NDT 
methods.  Use  of  the  finite  element  code  will  aid  in  the  development  of 
imaging  inverse  algorithms. 
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DIFFUSION,  WAVES.  PHASE  AND  EDDY 


CURRENT  IMAGING 

L.  Udpa  and  W.  Lord 

Electrical  Engineering  Department 
Colorado  State  University 
Fort  Collins.  CO  80523 

INTRODUCTION 

The  discovery  of  electromagnetic  induction  by  Faraday  and  Henry  In 
1831  not  only  served  as  the  catalyst  needed  for  the  very  creation  of 
electrical  engineering  but  also  provided  the  physical  basis  for  eddy 
current  nondestructive  testing  (NDT)  as  we  know  It  today  and  as  first 
realized  In  the  classical  experiments  of  Hughes^  ,  As  this  fundamental 
work  preceded  Maxwell's  prediction  of  electromagnetic  wave  phenomena  by 
over  half  a  century.  It  may  seem  somewhat  surprising  to  the  casual 
reader  that  there  should  be  any  need  to  explain  why  eddy  current  NDT 
phenomena  can  be  classified  as  quasl-statlc  in  nature  with  none  of  the 
attributes  of  classiced  electromagnetic  waves.  Unfortunately,  there  are 
many  misconceptions  concerning  the  wave-like  nature  of  eddy  current  NDT 
phenomena  which  have  even  led  to  the  suggestion^  that  conventional  eddy 
current  NDT  probe  signals  can  be  treated  holographically.  There  are 
several  reasons  for  the  existence  of  these  misconceptions: 

3 

1.  Many  papers  in  the  field  (see  for  example  Hochschild  ) 
describe  the  propagation  of  an  electronagnetic  plane  wave  in  a  medium  as 
being  analogous  to  eddy  current  NDT  phenomena.  Although  the  anadog 
itself  has  some  limited  validity,  it  is  rarely  if  ever  mentioned  that  a 
conventional  eddy  current  NDT  probe  does  not  launch  an  electromagnetic 
wave  (as  does  say  an  antenna) . 

2.  Solution  of  the  quasi-static  skin  effect  equation  for  current 
density  does  have  the  same  form  as  would  a  damped  electromagnetic  wave. 
However,  this  is  more  a  statement  of  the  consistency  of  Maxwell's  equa¬ 
tions  across  different  regimes  (see  Figure  1)  than  support  for  eddy 
current  waves.  A  number  of  authors  address  this  seemingly  anomalous 
situation  (see  for  example,  Stoll^  ,  Ferrari^  .  and  Melcher”  )  and 
clearly  differentiate  between  electromagnetic  diffusion  and  electromag¬ 
netic  wave  phenomena. 

3.  Much  of  the  terminology  associated  with  eddy  current  NDT 
phenomena  (phase,  for  example)  has  a  direct  counterpart  in  electromag¬ 
netic  wave  parlance. 
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Figure  1.  Overview  of  Maxwell’s  equation  in  different  frequency  regimes 


Conventional  eddy  current  NDT  phenomena  are  truly  ateady-state 
alternating-current  Induction  phenomena  completely  describable  by 
quasl-statlc  field  theory.  The  following  sections  of  this  paper  expand 
on  these  comments  and  discuss  their  implications  for  eddy  current  Imag¬ 
ing. 

PRINCIPLES  OF  EDDY  CURRENT  TESTING 

As  stated  In  the  introduction,  the  eddy  current  method  of  nondes¬ 
tructive  testing  Is  principally  based  on  Faraday's  law  of  electromag¬ 
netic  induction.  When  a  coll  excited  by  an  alternating  current  source 
is  brought  close  to  a  conducting  material,  the  primary  field  set  up  by 
the  coll  Induces  eddy  currents  In  the  material,  setting  up  an  opposing, 
secondary  field.  In  a  nonmagnetic  test  object,  this  results  In  a  reduc¬ 
tion  of  the  net  flux  linkages  of  the  coll,  thereby  reducing  the  Induc¬ 
tance  of  the  coll.  The  resistance  measured  at  the  terminals  of  the  coll 
Is  also  altered  to  account  for  the  eddy  current  losses  within  the 
material.  The  presence  of  a  defect  or  Inhomogenelty  In  the  material 
causes  a  redistribution  of  the  eddy  currents,  thereby  changing  the  com¬ 
plex  Impedance  of  the  probe  coll.  Changes  In  the  coll  Impedance  caused 
by  defects  In  the  material  are  represented  as  trajectories  In  the 
Impedance  plcUie  and  used  for  defect  characterization. 

From  considerations  of  the  operating  frequencies  and  '".mansions  of 
the  experimental  set-up.  the  eddy  currents  constitute  a  quasl-statlc 
phenomenon.  Under  these  conditions  the  displacement  current  Is 


neglected  and  Maxwell's  equations  are 

Vxi.-|f  (1) 

Vx  H  =  J  (2) 

V  •  B  =  0  (3) 

V  •  D  =  0  (4) 

Assuming  a  linear.  Isotropic  and  homogeneous  medium  the  constitu¬ 
tive  relations  are 

i  =  (1  H  (5) 

D  =  e  E  (6) 

J  =  o  i  (7) 


Decoupling  equations  (1)  and  (2)  using  the  constitutive  relations,  the 
governing  equations  for  the  fields  and  currents  are 


\r  j  -  at 


(10) 


For  single  frequency  sinusoidal  excitation,  the  one-dlmenslonal 
form  of  (10)  can  be  written  In  the  phasor  form  ais 


=  jupeJ 
dx 

(11) 

The  steady,  state  solution  for  the  case  of  a  sheet  of  alternating 
current  density  at  the  surface  of  a  seml-lnflnlte  medium  is  given  by 

J(x,t)  =  J(0,t)exp(-|)exp(-J(|  -  »t)} 

(12) 

where  the  skin  depth  6  is  given  by 

(ittlO 

(13) 

Though  equation  (12)  has  the  mathematical  properties  of  an  attenuated 
wave  it  does  not  describe  a  true  physical  wave.  What  it  actually 
describes  is  the  distribution  of  steady  state  alternating  currents  whose 
magnitude  decays  exponentially  and  whose  phase  atngle  varies  linearly 
with  depth  according  to  the  relations 

J(x)  =  J(0) 

(14) 

9(x)  = 

(15) 

Equations  (12)  -  (15)  are  derived  for  a  rather  contrived  geometry  for 
the  sake  of  computational  ease.  The  eddy  current  paths  In  an  actual 
eddy  current  test  au'e  far  more  complex  and  the  presence  of  anomalies  In 
the  medium  further  perturb  this  distribution.  However  the  solutions 
obtained  above  do  serve  to  Illustrate  the  general  behavior  of  fields  and 
currents  and  the  associated  skin  depth  which  is  of  significant  impor¬ 
tance  in  eddy  current  NDT  situations. 

The  phase  and  magnitude  described  above  for  a  prescribed  depth  are 
not  experimentally  measurable.  VThat  is  measured  In  an  eddy  current  test 
Is  the  complex  Impedance  of  the  probe  coll,  which  is  affected  by  the 
total  eddy  current  distribution  in  the  test  specimen. 

The  measured  Impedance  of  the  coll  can  be  expressed  as 

^0  ■  "o  *  (I«, 

where  V  and  I  are  RMS  valued  voltage  and  current  phasor s. 


The  phase  angle  of  the  eddy  current  test  data  is  then  given  by 

-1  *c 
d  -  Tan  ^(^) 

°  ”c 


which  does  not  have  a  simple,  direct  relationship  with  the  6  of  equation 
(15). 


The  eddy  current  test  can  be  compared  to  the  operation  of  a 
transformer  where  the  probe  coil  is  the  primary  and  the  test  material 
constitutes  the  secondary.  Just  as  in  a  transformer,  the  secondary  pro¬ 
perties  are  referred  to  the  primary  side  and  hence  the  material  charac¬ 
teristics  are  reflected  in  the  coll  impedance  measurement.  However  this 
measurement  reveals  nothing  about  the  actual  current  distribution  within 
the  test  specimen.  Determination  of  materiail  characterstics  on  the 
basis  of  the  test  signal  from  the  eddy  current  probe  is  therefore  a  com¬ 
plex  process. 

IMAGING 

Imaging  or  inversion  of  eddy  current  data  is  the  problem  of  recon¬ 
struction  of  the  defect  in  three  dimensions,  given  the  meaisured  signal. 

A  direct  approach  by  use  of  a  theoretical  model  based  on  the  underlying 
physical  process  is  too  complex  and  most  of  the  existing  defect  charac¬ 
terization  schemes  resort  to  the  indirect  algorithmic  methods  which 
depend  on  characteristic  features  in  the  signal  for  classification 
information. 

I 

HOLOGRAPHIC  IMAGING 

In  an  attempt  to  directly  image  the  defect  in  three  dimensions, 
Hildebrand  et  al.^  apply  holographic  principles  to  eddy  current  data, 
interpreting  the  eddy  current  phenomenon  as  an  interference  between 
incident  and  reflected  electromagnetic  waves.  The  magnitude  and  phase 
of  the  coil  impedance  data  in  (16)  are  thus  Interpreted  as  the  magnitude 
and  phase  of  a  scattered  wavefront  that  satisfies  the  Helmholtz  wave 
equation.  The  method  then  applies  a  backward  wave  propagation^  algo¬ 
rithm  to  the  eddy  current  data  to  reconstruct  the  defect  in  three  dimen¬ 
sions  . 

This  procedure  gives  meaningful  results  only  if  the  data  input  to 
it  indeed  describes  a  true  wave.  Otherwise,  the  method  functions  as  a 
low  pass  filter  with  a  phase  response  given  by 

{l-(Xu)^-(l,v)^}^^^2]  which  merely  distorts  the  eddy  current 
probe  ^gnals. 

ALGORITHMIC  IMAGING 

Algorithmic  imaging  is  a  procedure  for  deducing  defect  geometry 
parameters  by  using  distinctive  properties  of  the  measured  signal.  The 
algorithmic  methods  for  characterizing  and  sizing  defects  include  signal 
processing  techniques  such  cis  adaptive  learning  networks^  ,  and  the  use 
of  Fourier  descriptors  developed  by  Udpa  and  Lord®  .  In  edl  these 
methods,  the  eddy  current  probe  signals  are  treated  as  signatures  of  the 
defect  that  produced  them.  The  signal  from  each  defect  type  is  t^en 
represented  by  a  set  of  features  either  from  the  time  domain  or  fre¬ 
quency  domain  or  a  combination  of  both.  A  data  bank  of  feature  vectors 
corresponding  to  all  the  expected  defect  types  is  thus  built.  An  unk¬ 
nown  signal  is  then  classified  as  belonging  to  one  of  the  sets  in  the 
data  base  by  pattern  recognition  techniques. 


A  more  direct  approach  has  been  used  by  Copley^  where  the  phase  and 
amplitude  of  the  measured  signal  are  directly  interpreted  in  terms  of 
defect  dimensions  by  use  of  a  set  of  calibration  curves. 

For  the  sake  of  comparing  the  results  of  holographic  imaging  and 
algorithmic  imaging  a  simple  image  processing  algorithm  was  implemented. 
Using  the  experimental  set  up  shown  in  Figure  2.  the  horizontal  and 
vertical  channels  of  the  complex  probe  data  were  sampled  at  discrete 
spatial  points  digitized  and  stored  on  the  VAX  11/780  computer.  The 
aluminum  bar  containing  defect  patterns  such  as  that  shown  in  Figure  3 
was  used  as  the  test  specimen.  Scanning  was  done  on  the  other  side  of 
bar.  so  that  these  90%  through-wall  holes  served  as  subsurface  defects. 

In  Figure  4  the  complex  data  is  displayed  as  four  different  grey-level 
coded  images  representing  vertical  channel  data,  horizontal  channel 
data,  magnitude  and  phase.  The  result  of  holographic  imaging  is  shown 
in  Figure  S  where  the  probe  data  is  ’back  propagated’  to  the  plane  of 
the  defect.  The  result  of  a  basic  thresholding  and  edge  detection  algo¬ 
rithm  is  shown  in  Figure  6.  This  algorithm  also  computes  the  diameters 
of  the  holes.  Space  limitations  preclude  a  complete  discussion  of  adl 
the  test  results  comparing  holographic  and  algorithmic  imaging  results. 
The  reader  is  left  to  draw  his  own  conclusions  from  Figures  5  and  6. 

CONCLUSIONS 

Holographic  imaging  algorithms  can  certainly  be  applied  to  the 
analysis  of  conventional  eddy  current  probe  data.  The  key  questions 
highlighted  by  this  paper  are: 

1.  Does  such  a  procedure  make  any  physical  sense? 

2.  Does  such  a  procedure  have  distinct  advantages  over  algorithmic 
imaging? 

In  answer,  eddy  current  phenomena  are  quasi-static  phenomena  where 
the  operating  frequencies  and  characteristic  dimensions  are  such  that 
the  displacement  current. is  negligible.  Consequently  the  eddy  currents 
are  described  by  diffusing  phasors  rather  than  by  attenuated  waves. 
Secondly,  in  an  eddy  current  test  the  measured  probe  impedance  is  caused 
by  the  integrated  effect  of  adl  the  currents  in  the  specimen  and  hence 
the  phase  of  the  test  data  cannot  represent  the  phaise  of  induced  eddy 
currents  or  the  phase  of  an  electromagnetic  wave.  Thus  the  indirect 
nature  of  the  eddy  current  test  makes  a  direct  approach  to  the  inverse 
problem  very  complex  and  algorithmic  methods  appear  to  be  a  more  practi¬ 
cal  solution. 
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Figure  3.  A1  bar  with  defect  patterns. 
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INTRODUCTION 

Numerical  analysis  techniques  have  been  successfully  applied  to  the 
modeling  of  electromagnetic  field/defect  Interactions^  .  Studies  of 
magnetostatic  leakage  field  and  eddy  current  NDT  phenomena  have  clearly 
shown  that  finite  element  codes  can  be  used  effectively  for  probe 
design^  and  the  simulation  of  test  geometries  difficult  to  replicate  In 
the  laboratory^.  In  extending  these  codes  to  three  dimensional 
geometries^  and  pulsed  eddy  current  phenomena  ,  it  was  realized  that  the 
required  computing  capability  should  also  be  sufficient  to  model 
ultrasound/defect  interactions  directly  in  the  time  domain.  Increasing 
availability  of  powerful  vector  computers”  bodes  well  for  the  ultimate 
solution  of  the  generic  NDT  problem  in  which  it  is  desired  to  predict 
the  probe  response  to  any  arbitrarily  shaped  defect.  As  a  first  step  in 
this  direction,  the  NDT  research  group  at  Colorado  State  University, 
following  the  pioneering  numerical  efforts  of  Bond'  and  Dewey”'  has 
developed  a  finite  element  code  for  direct  time  domain  solution  of  the 
elastic  wave  equatipn  (Figure  1  shows  the  relationship  between  numerical 
auid  analytical  approaches) .  The  following  sections  describe  the  finite 
element  formulation  and  the  application  of  the  code  to  the  prediction  of 
2-D  displacements  In  a  rectangular  bar  excited  at  one  end  by  a  step 
input  of  force. 

FORMULATION 


The  general  equation  of  motion  can  be  written  in  the  form 

7.T+F-p^ 

at^ 


(1) 


where  T.F.u  represent  stress  tensor,  body  force  and  displacement  vectors 
respectively,  p  denotes  the  material  density.  Three  restrictions  are 
imposed : 

1)  no  body  forces 


F  -  0 


2)  no  Internal  energy  losses  and  small  deformations  such  that 
Hook’s  law  Is  applicable 


T  »  C:S 

with  C  being  the  forth  rank  material  tensor  and  S  representing 
the  strain  tensor 

3)  only  a  homogeneous  Isotropic  solid  Is  considered.  Thus,  the 
material  tensor  consists  of  only  two  Independent  coefficients  X 
and  4  (Lame  constants) 

^Ijki  “  ^It^jk^ 

Substitution  of  these  three  conditions  Into  (1)  yields  the  elastic  wave 
equation  In  rectangular  coordinates 

2- 

(li+X)v(7«u)  +  pV^u  =  p  ... 

dt^ 

2  2 

If  V,  s  (X+24)/p  and  V  =  p/p  are  Introduced  as  longltudlnad  and  shear 
velocities,  (2)  can  be^expressed  for  the  two  dimensional  case  as 

2  a^u  d\  d\  d\ 

^  3x2  ^s  3y2  axay  ^^2  (3a) 


v2  2  aa!]^  ^  2_^  ^ 

^  ay2  «  ax2  L  ayax  ^^2 


with  the  Neumann  type  boundary  conditions  given  by 

^  .  au 
T  •»  ^  +  (v2-2v2)-— ^ 
xx/p  V  ax  ay 


-  au  au„ 

r  *  v2( — ^  +  — ^)  =  T 
xy/p  s' ay  ax'  V/p 


yy/p  L  ay 


«  ..  ..  3u 

vf  -rf  +  (V?-2v2)-^ 
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FINITE  ELEMENT  IMPLEMENTATION 


Instead  of  developing  a  direct  discretization  of  (3)  by  means  of 
collocation  or  Galerkln's  method,  we  consider  an  energy  related  func¬ 
tional 


z  L  dx  oy 

2  2  2  ®“v  2 

*  ^  ♦  -if*  «<“,  -f*",  rf 

OV  «w 


)}p<iv 


which,  upon  finding  a  stationary  value  with  respect  to  the  unknown  dis¬ 
placements  u  .u  ,  results  in  the  same  solution.  An  easy  way  to  check 
the  correctness^of  the  above  functional  Is  to  utilize  variational  cal¬ 
culus  in  order  to  cu:*rlve  at  the  so  called  Euler  equations  which  subse¬ 
quently  yield  the  original  elastic  wave  equations  (3a)  and  (3b).  It  can 
also  be  shown  by  the  same  derivation  that  the  stress  free  boundary  con¬ 
ditions  are  implicit  in  the  energy  related  functional. 

To  solve  (S)  in  terms  of  the  unknown  displacements,  the  following 
four  steps  have  to  be  performed 

a)  discretize  solution  domain  into  a  finite  number  of  elements 

b)  find  a  stationary  value  for  (S)  with  respect  to  u  .u 

X  y 

c)  replace  u  .u  by  the  approximations 
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J  etc.,  where  [N(x.y)l 
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denotes  the  shape  functions  as  a  row  vector  with  (u  )  .  (u  }  being 
the  unknown  displacements  at  the  nodal  points  of  eaSb^elemeXt?  The 
resulting  elemental  matrix  equation  takes  on  the  form 
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with  the  coefficients  of  the  submatrices  given  by 
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F  (I) .F  (I)  are  external  driving  forces.  The  numerical  Integration 
Is  carried  out  by  employing  a  7  point  Gaussian  quadrature  formula. 


d)  Assemble  all  the  elemental  matrices  (6)  Into  a  global  matrix  which 
can  be  solved  for  (u  }  and  {u  )  Before  the  assembly  can  be  done 
however,  the  problem^conslsts  of^ Integrating  the  second  time  deriva¬ 
tive  In  (6).  Possible  Integration  schemes  are  tez*med  as  either 
explicit  or  Implicit  depending  on  whether  a  matrix  Inversion  of  [Kl 
Is  Involved.  The  central  difference  Integration  (explicit)  as  well 
as  the  Houbolt.  Wilson  and  Ne%imark  Integration  (Implicit)  have  been 
Implemented.  For  the  purposes  of  this  paper  only  the  Newmark 
Integration  Is  given 
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This  scheme  can  be  made  unconditionally  stable  depending  on  the 
selection  of  a  and  6. 


APPLICATIONS 


In  order  to  validate  the  finite  element  code,  a  bar  subject  to  a 
step  tension  T.  was  modeled  as  shown  In  Figure  2.  In  Figure  3  the  u 
displacement  Is  plotted  at  three  different  locations  A.B.C  within  the 
bar.  The  results  are  In  excellent  agreement  with  the  one  dimensional 
displacement  predictions  by  Dewey  et  al.  who  also  shows  the  analytical 
series  solution. 


A  more  crucial  test  results  If  one  attempts  to  obtain  the  u  dis¬ 
placements  for  a  wide  rectangular  bar  subject  to  a  longitudinal  Step 
pressure  loading.  Jones  and  Ellls^  compared  the  theoretical  predictions 
of  the  plane-stress  theory  with  their  experimental  observations.  Their 
experimental  results  for  a  130  Inch  long  and  1.5  Inch  wide  bar  with 
V  =1.248  X  10^  In/s  and  a  Poisson  ratio  of  O-0.335  also  show  good  agree- 
ffllnt  with  the  finite  element  prediction  shown  In  Figure  4.  Typical 
solution  times  on  a  VAX  11/780  computer  for  two  different  mesh  sizes  are 
given  In  Table  1.  These  figures  are  somewhat  misleading,  however,  since 
both  mesh  size  and  computer  code  €U'e  not  yet  optimized.  In  general,  the 
computer  time  Is  a  function  of  wave  velocity,  transducer  frequency,  sam¬ 
pling  rate,  distance  of  travel  and  mesh  size.  To  Illustrate  this  for  a 
two  dimensional  case  and  show  how  a  powerful  computer  like  the  CYBER  205 
can  significantly  reduce  the  time  requirement,  consider  the  more  sophis¬ 
ticated  problem  of  pulse  echo  wave  propagation  as  shown  In  Figure  5 . 

Here  a  pulse  created  by  a  1  MHz  transducer  propagates  with  a  longitudi¬ 
nal  velocity  of  V  =5000  m/s  through  a  specimen  of  12.5  cm  thickness. 

The  total  travel  time  for  twice  the  thickness  Is,  therefore.  50)is. 

Based  on  an  assumed  sampling  rate  of  32  MHz  It  follows  that  1600  time 
steps  solution  time  are  required.  If  a  solution  domain  requiring  a  mesh 
size  of  3000  nodes  Is  assumed.  It  will  take  a  VAX  11/780  computer  about 
4.5  minutes  to  solve  the  resulting  matrix  equation  at  each  time  step. 
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Fig.  1.  Relationship  between  nuaier-  Fig.  2.  Geometry  and  boundary  condl 
leal  (finite  differences,  tlons  of  rectangular  bar 

finite  elements)  and  analyt¬ 
ical  approaches 


g.  3.  Normalized  u^  displacement  Fig.  4.  Normalized  u^  displacement 
at  locations  A,B,C  at  x-37.3  In. 
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CP— requirements  on  VAX  1 1  /78Q 
elements  I  nodes  |  CP-time(900  timesteps) 


Based  on  our  experience  with  eddy  current  calculations,  the  CYBER  205 
reduces  the  solution  time  to  O.OOSd  min.  For  a  total  of  1600  time  steps 
this  would  amount  to  120  hours  on  the  VAX  versus  9  minutes  on  the  CYBER 
205. 

CONCLUSIONS 

A  considerable  amount  of  work  remains  to  be  done  before  numerical 
code  can  be  used  as  an  engineering  tool  for  the  design  amd  anedysls  of 
ultrasonic  nondestructive  tests.  Early  studies  in  this  field  show  prom¬ 
ise.  however,  and  the  increasing  availability  of  supercomputers  can  pro¬ 
vide  the  computational  power  needed  to  ultimately  predict  ultrasonic 
transducer  responses  from  realistic  defect  geometries.  Although  the 
finite  element  formulation  and  applications  described  in  this  paper  are 
2-D  in  nature,  only  computational  cost  limits  the  extension  to  3-D 
geometries . 
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